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Abstract. The morphological distribution of primordial den¬ 
sity peaks is assessed. Previous determinations, those of Pea¬ 
cock & Heavens (1985) and those of Bardeen et al. (1986) have 
contradictory concluded that there exist a tendency towards 
prolate or towards oblate shapes, respectively. By using two 
methods, the Hessian and the inertia tensor momenta we have 
performed numerical determinations of the triaxiality of den¬ 
sity perturbations in Gaussian random fields with power law 
spectra. We show that there is no present any tendency of 
shapes, and that the triaxiality distribution is independent of 
the spectral index. Moreover, it is shown that the results of 
Peacock and Heavens are compatible with our determinations. 
These results are in complete agreement with current triax- 
ial distributions infered for galaxies and clusters both, from 
observations and from numerical simulations. 
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1. Introduction 

In the extensive analysis of Random Gaussian fields by 
Bardeen, Bond, Kaiser & Szalay (1986, BBKS), it was found 
that the density maxima are intrinsically triaxial with some 
preference towards oblate shapes. Several authors have at¬ 
tempted to connect such a distribution of triaxiality, or the 
height of the peaks in the primordial density field to the ob¬ 
served distribution of triaxiality of galaxies and clusters. For 
instance, Dubinski (1992) based on the Peacock & Heavens 
(1985, PH) results, compared the triaxiality distribution of 
density peaks to the triaxiality of dark halos in hierarchical 
models. He aimed to study whether the initial triaxiality dis¬ 
tribution, with a tendency towards prolateness, is preserved 
during the non-linear evolution of the density field. His hypoth¬ 
esis was, that if the shape of peaks indeed determines the shape 
of dark haloes around galaxies, an excess of prolate dark halos 
may be expected because of the natural bias for prolate initial 
conditions, which he did not find. Furthermore, Evrard (1989) 
and Evrard, Silk & Szalay (1990) tried to explain the Hub¬ 
ble sequence by identifying elliptical galaxies with the highest 
peaks in the density field, and spirals with the lowest peaks. 
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Some other authors have attempted to infer the 3- 
dimensional shape of galaxies and clusters from the observed 
projected axial ratios. Fasano & Vio (1991), for example, anal¬ 
ysed 204 elliptical galaxies and showed that the distribution of 
triaxiality of these galaxies is incompatible with both, a dis¬ 
tribution of pure oblate or pure prolate spheroids. A very low 
number of round elliptical galaxies was found in their sample. 
Lambas, Maddox & Loveday (1992) studied a larger sample of 
20399 galaxies, including spiral and ellipticals galaxies. They 
confirmed that elliptical galaxies are triaxial with no preference 
for oblate or prolate shapes. 

Clusters of galaxies have also been studied by e.g. Carter & 
Metcalfe (1980), who found not only that clusters of galaxies 
are flatter than elliptical galaxies, but also that their distribu¬ 
tion of shapes is better fitted with a set of oblate spheroids. On 
the other hand, Binggeli (1982) analysed the apparent elliptic- 
ity of 44 Abell clusters, finding that the shapes distribution is 
consistent with prolate shapes. Furthermore, Plionis, Barrow 
& Frenk (1991), with a more complete sample of 6000 clusters 
determined that the morphological distribution is inconsistent 
with pure oblate objects, and conclude that a triaxial distri¬ 
bution of shapes would fit better their sample. Such a trend 
has been confirmed by using numerical simulations in hier¬ 
archical models by Efstathiou et al. (1985) and Frenk et al. 
(1988). Concerning superclusters, there also exist contradic¬ 
tory claims regarding their shapes. West (1989) generated a 
catalog of 48 probable superclusters by using a large sample 
of Abell clusters and found that their shape is consistent with 
prolate objects. Meanwhile, Plionis et al. (1992) by employing 
several statistical techniques and a sample of Abell and ACO 
(Abell, Corwin & Olowin, 1989) clusters, obtained that super¬ 
clusters are triaxial with a strong preference toward oblatness. 
Because the density field is in the linear regimen on large scales, 
it is possible that the morphology of large-scale perturbations 
shall reflect the primordial shapes. In adition, there are theo¬ 
retical evidences in favor that information of the initial shape 
at the moment of galaxy formation can be preserved after the 
relaxation process (Aarseth & Binney, 1978). 

Clearly, there is no consensus on the distribution of tri¬ 
axiality for the different observed population of cosmic ob¬ 
jects. Whilst several statistical properties of Gaussian and Non- 
Gaussian statistics have been addressed by e.g. BBKS, Couch- 
mann (1987), Barrow & Coles (1987, 1990) and Coles (1990), 
concerning shapes of peaks in Gaussian fields there are only two 
studies with contradictory conclusions. For PH peaks present a 
tendency towards prolate shapes, whereas the analytical study 
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by BBKS gave a trend towards oblate shapes. The real shape 
and triaxial distribution of density perturbations, are highly 
desirables if one aims to understand the effect of the non-linear 
evolution of the density field on the morphology of perturba¬ 
tions, and moreover, the effect of tidal interactions. In this 
paper, we perform numerical determinations of the triaxiality 
of density peaks and show that Gaussian random fields present 
no tendency for oblate or prolate shapes. They are intrinsically 
triaxials. We also confirm the BBKS statement that triaxial¬ 
ity distribution is not sensitive to the spectral index of power 
law spectra and therefore is useless to constraint cosmological 
models. For these goals the rest of the paper is organized as 
follows. In Section 2, we start by briefly review the theretical 
procedure to assess shapes. A description of the adopted nu¬ 
merical methods is done in Section 3. Section 4, contains our 
results and a comparison with prior determinations. A discus¬ 
sion and our conclusions are presented in section 5. 


P(e,p | x) = (^y) 1 / 2 J^ exp(-(5/2)x 2 (3e 2 +p 2 )))W(e,p),(6) 

where we have used x = x* = 717 the most probable value 
of the curvature, and 7 is a function of the spectral index 
7 = [(n + 3 )/(n + 5)] 1/<2 . Hereafter the height of the peaks 
will be expressed in units of the root mean square <r of the 
density field as <5 = ua. Thus, given a probability P and a 
threshold lieigth v max of the peaks, equation ( 6 ) defines an 
isoprobability curve which encloses at least 100(1 — P)% of all 
the peaks in the density field with height v < v max . Figure 
1 exhibits some of the isoprobability contours as a function 
of the probability and height of the peaks for the spectral in¬ 
dexes n=l, —2. Only these spectral indexes were included to 
enhance any possible dependence of the triaxiality distribution 
on the spectrum index. Nevertheless, no extreme dependence 
is observed. Our numerical results will be compared with the 
theoretical contours of Figure 1. 


2. Shape of Peaks: theoretical formulation 

In the following analysis we shall adopt the procedure and 
notation of BBKS. To first approximation one can expand the 
density field around a peak as 


5(r) = <5(0) + r 


d5(r) 1 d 2 S(r) 

~drT r=0 + 2 nri 'd^d7 ] 


r=o +..., 


(1) 


where the second term on the r.li.s. is zero due to the ex¬ 
tremum condition. The term composed by the second deriva¬ 
tives of the field is the Hessian matrix 


H=- 


d 2 S 

dridrj ’ 


( 2 ) 


whose eigenvalues ordered, without lack of generality, as 
Ai > A2 > A3 represent the quadratic approximation 


<5(r) = 5(0) ~y^Ajr 2 /2, 


(3) 


which defines an ellipsoid with semiaxes 


_ r 25(0) — 5t ll /2 
“ 1 Ai J ’ 


(4) 


3. Numerical Approximation 

We have used the Hessian and the inertia tensor methods to 
assess the intrinsic shape of peaks. THe first one, was used by 
PH to carry out numerical determinations of the morphology 
of peaks, whereas the second method was adopted by BBKS 
for their analytical treatment. Since we aim to numerically test 
and collate these results, both methods require of smoothing 
the discrete density field, randomly generated in a 64 3 box, and 
an interpolation process. We based this latter on the triangular¬ 
shaped cloud method (Hockney & Eastwood, 1988). This en¬ 
ables us the reconstruction of the density field at any point 
within the cubic grid. Further, in order to apply the interpo¬ 
lation, it is also required to set up the minimum number of 
contributing points-density to the inertia momenta necessary 
to get stable results of the morphology. In the calculations, we 
have associated a unit mass to each point of the peak, though 
various others weighting schemes can be used, see e.g. West 
(1989b) and Plionis et al. (1992). A fully detailed description 
of the method used, and the determination of the minimum 
number of grid points and steps of interpolation can be found 
in Gonzalez (1994). N s denotes the number of sub-divisions of 
each l 3 cube that is crossed by the isodensity surface. 


where 5t defines the value of the isodensity countour for 
which the triaxiality is to be determined. Equation (4) is the 
definition adopted by PH to analyse the morphology of density 
peaks in terms of the ratio of their main axes. 

Aditionally, the triaxiality parameters are defined by BBKS 
in terms of the Hessian eigenvalues as 

Ai — A3 Ai — 2 A 2 + A 3 

2EV ’ (5) 

where e measures the ellipticity in the Ai — A 3 plane and p 
measures the prolatness or oblatness. If 0 > p > — t then the 
ellipsoid is prolate-like, while if t > p > 0 it is oblate-like. The 
limiting cases are, p = —t for prolate spheroids and p = e for 
an oblate spheroid. 

The conditional distribution, P(t,p \ x), which describes 
the probability of finding a density perturbation with triaxial¬ 
ity parameters e and p in the interval e + de and p + dp, given 
the value of the normalized curvature x = V 2 5 /ct 2 (see eq. [7.6] 
in BBKS), is 


3.1. Methodology 

We generated random fields with a power law spectrum, 
P(k) = Ak n , with spectral indexes n = 1,0, —1,-2. The 64 3 
box used in the calculation has a physical length of 64h~ 1 Mpc 
on a side. We then proceed as follows: 

1. The density field is convolved with a Gaussian filter 

W G = exp(—(k 2 R 2 )/2), (7) 

where k is the perturbation wavenumber and R is the fil¬ 
tering radius. 

2. The positions of grid maxima are found, which are defined 
as those grid points with overdensity S max higher than 
the density of their 26 nearest neighbours. In general, the 
grid maxima positions do not coincide with position of real 
peaks. 

3. The value of the isodensity surface 5t, for which the triaxi¬ 
ality parameters will be calculated, is chosen by introducing 
a factor / 

St = fSmax, ( 8 ) 
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Fig. 1. Isoprobability contours of the distribution of triaxiality for peaks as a function of their height is, for n = 1, —2. The theoretical 50% 
isoprobability contours are included, which shall contain half of the total number of peaks. 


with 0 << / < 1. Then, the inertia tensor is calculated 
1 m 

hi = — yTfffc- < Xk >)(xi- < XI >)], (9) 

i= 1 

which involves the use of the grid points and the interpola¬ 
tion code of N s = 14 steps, and where k,l = 1, 2, 3, Xu, X 2 i 
and X 3 i are the Cartesian coordiantes of the ith point. For a 
given /, the sum runs over all the cubic nodes with 8 < St- 
The inertia tensor eigenvalues a§ > a| > a?, are then re¬ 
lated to the eigenvalues of the Hessian matrix 

A, = 1/a?, A 2 = 1/a? and A 3 = 1/a?. (10) 


4. In order to assure the validity of the quadratic approxima¬ 
tion, the calculation of the Hessian matrix must be made 
close to the position of the real peak. Therefore, one of 
our requeriments of the inertia tensor method, is to per¬ 
form an interpolation process to localize the real maxi¬ 
mum. Afterwards, the Hessian matrix is constructed by 
interpolating the density at points separated at distances 
d = R/ 3, R/2, R from the peak center. These three dis¬ 
tances were used to test the numerical stability and con¬ 
stancy of the shape around peaks. 
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Fig. 2. Morphology of peaks as a function of their heigh, for n = 1. The theoretical 50% isoprobability contours is included, which contains 
approximately half of the total number of peaks. The mean triaxiality parameters are < e >= 0.174 and < p >= 2.3 X 1 0" 2 . 


4. Results 

4-1. Statistical analysis 

Figure 2 shows the distribution of the triaxiality parameters as 
a function of the threshold height. It confirms the trend for the 
highest peaks, v > 3.5, to be more spherical than the lower ones 
as found by BBKS. The large dispersion of peaks shapes for all 
over the e — p plane, comes from the peaks lower than v < 1.5cr. 
Meanwhile, the distribution of shapes for peaks in the interval 
1.5 < v < 3, -where the number density of peaks is maximum- 
, is more concentrated towards sphericity and at least 50% of 
them lay within the theoretical isoprobability contours. If we 
look at the form of these isoprobability contours, we will notice 
that the major part of the area bounded by them, is into the 
oblate shapes side. In order to check whether this should be, 
or not, understood as a tendency towards oblateness, we have 
calculated the ratio F of the number of oblate ellipsoids N 0 to 
the number of prolate ones N p , both of which are numerically 
obtained. Any significant trend must be reflected in F. In the 
second and third columns of Table 1 we display this ratio for 
different spectral indexes as a function of the height of peaks. 
It is easily oberserved that there is not any morphological pref¬ 
erence. Moreover, when we only consider the highest peaks no 
strong trend is detected, even when they constitute a small 
sample. 

The histogram of Figure 3 displays the distribution of 
oblate and prolate ellipsoids. This histogram and the ratios 
of Table 1, suggest that if an especific trend of shapes exists, 
in fact it behaves stochastically, and has no statistically rel¬ 
evance. Therefore, an intrinsic triaxial shape of density per¬ 
turbations would be the most acceptable conclusion. An in¬ 
spection of the distributions of Figure 3 also confirms that the 
results are not significantly dependent on the spectral index. 
The distribution of ellipticity and prolateness for all the models 
exhibit a slight tendency for oblate shapes. The mean values 


< t >= 0.18 ± 0.02 and < p >= (.13 ± 0.003) x 10 -3 incorpo¬ 
rate all the models. These values and the distribution of peaks 
in the e — p plane show disagreement both with the theoreti¬ 
cal predictions of BBKS and with the numerical ones of PH. 
In order to compare our results with those of these latter au¬ 
thors, it is necesary to assess the shapes by resorting to the 
Hessian matrix as indicated in section 3. Once the Hessian is 
diagonalized we take the ratio of axes defined by equation ( 10 ) 


/AsA 1 / 2 _ ai _ short — axes 

V Ai / 03 long — axes 

/ AsX 1 / 2 02 middle — axes 

V A 2 / 03 long — axes 


( 11 ) 

( 12 ) 


Figure 4 is an example of the distribution of the ratio of axis 
in the s — m plane. A visual judgement of such a distribution 
can easily lead us to wrongly conclude -as PH did that a 
tendency of the peaks exists towards prolateness. In a s — rn 
diagram, a difference between a slight tendency to oblatness 
and a slight tendency to prolateness would be quite clear and 
not neglegible. The reason for this artificial discrepancy with 
our results is that PH did not consider the line which separates 
these two tendencies: p = 0, V e in the interval 0 < e < 0.33, 
i.e. all the eigenvalues which satisfy 


Ai — 2 A 2 + A 3 = 0, 


(13) 


or in terms of the ratio of axes 



which gives the equation of the curve 
m 

y/1 — m 2 


(14) 


(15) 
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Shapes of Density peaks 

Spectral index 

l st Realization 

2 n<1 Realization 

3 m Realization 

n 

(Inertia Tensor) 

(Inertia Tensor) 

(Hessian) 


r = N 0 /N p 

r = n 0 /n p 

r = No/N p 

n = 1 

v > 3.5 

110/105 

112/111 

104/101 

v > 2.5 

499/495 

502/506 

503/498 

v > 1 

3140/3144 

3093/3092 

3098/3101 

v > 0 

3283/3276 

3337/3341 

3305/3306 

n — 0 
v > 3.5 

93/95 

113/109 

107/102 

v > 2.5 

413/409 

402/410 

423/417 

> 1 

2527/2512 

2314/2321 

2323/2321 

v > 0 

2687/2693 

2498/2475 

2548/2539 

n = —1 
v > 3.5 

70/68 

73/76 

69/72 

v > 2.5 

297/305 

305/308 

313/309 

v > 1 

1789/1797 

1810/1797 

1805/1812 

v > 0 

2057/2048 

2023/2026 

2063/257 

n = —2 
v > 3.5 

37/32 

42/41 

38/41 

v > 2.5 

153/158 

165/163 

157/155 

v > 1 

1071/1083 

1112/1096 

1086/1083 

v > 0 

1386/1394 

1367/1356 

1352/1343 


Table 1. Ratio of the number of oblate objects to the number 
of peaks. The used method is indicated. 


When this function is included in the s — m diagram of 
Figure (4), it is observed that the results of PH are also consis¬ 
tent with no tendency, perturbations are intrinsically triaxial 
as is further shown both by the distribution of the number of 
oblate and prolate shapes of Figure 5 and by the fourth col¬ 
umn of Table 1. Going back to the results of BBKS one would 
wonder why there should be a natural bias for oblate shapes ? 

4-2. Random generation of Ai 

The primordial density field is a Gaussian random field, i.e. it is 
generated by superposing a large number of Fourier modes with 
phase angles drawn at random from a uniform distribution in 
the interval 0 — 2n. The shapes of density peaks, regarding them 
since the Hessian point of view, are the result of the curvatures 
of the density field on the three main axes. The eigenvalues of 
the Hessian matrix are therefore randomly generated numbers 
as well. The value of one does not affect the value of the other 
two. Under this assumption, it is no clear how a tendency to¬ 
wards oblateness can arise. It is also under this assumption 
that we have generated triplets of random numbers, without 
paying attention for them to fulfill the basic requirements of a 
cosmological density field; power law spectra, spectral index, 
etc. This has the goal of trying to reproduce our triaxiality 
distribution as a random process. This would indicate that the 
tendency infered from the BBKS results is a consequence of 
several analytical simplifications which hide the real triaxial 
shape of peaks. 

We will call the triplets, eigenvalues,such that 0 < A* < 1. 
Once they are put in an increasing order, we determine the tri¬ 
axiality parameters through equation (5). The results of this 


of prolate one, as a function of the spectral index n and the high 


Fig. 3. Distribution of oblate and prolate shapes of peaks as a func¬ 
tion of the spectral indexes, n = 1,-2, and hight v > 0. These 
distributions correspond with the cases of Table 1. The distribution 
in dashed line is a theoretical distribution of triaxiality, constructed 
by generation of random triplets A, (see Sec. 4). 


are also displayed in the distribution of Figure 3 (in dashed 
line), which present no substantial differences with the distri¬ 
bution of shapes for a Gaussian random field obtained by using 
the inertia tensor and Hessian methods. Moreover, we can pro¬ 
vide a meassure of how equal these two distributions are, in 
terms of their moments; mean < p >, standard deviation a, 
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Fig. 5. a) The same as Figure 4 for n = — 1, but here the boundary curve (eq. 15) between the oblate and the prolate shapes is included, 
b) This distribution of oblate and prolate shapes, and fourth column of Table 1, show that the intrinsic shape of peaks is triaxial 


Fig. 4. Distribution of peaks according with the ratio of axis, s and 
m, as presented by Peacock and Heavens (1985), where an artificial 
tendency towards prolate shapes is observed. 


skewness s and kurtosis k. These moments characterize the ob¬ 
served asymmetry of the triaxiality distribution. For the Gaus¬ 
sian density field we obtain: p = 1.324 x 1CU 2 , cr = 0.1189, s = 
0.9091 and k = 0.9364. Whereas for the distribution deter¬ 
mined by the random triplets we get: < p >= 8.53 x 1CU 3 , a = 
0.1194, s = 0.9251 and k = 0.9635. Therefore, it is quite likely 
that this two distributions are generated by the same stochastic 
process. 


5. Conclusions 

The present investigation underlines the nature of primordial 
density peaks and the dependence of the triaxiality distribution 
on the spectral index, for power law spectra. This numerical 
analisys explored the veracity of the PH and BBKS conclusions 
on these points. 

For all the spectral indexes, we confirmed the BBKS re¬ 
sult that high peaks tend to be more spherical than the lower 
ones. When the whole sample of peaks is considered, an intrin¬ 
sic triaxial morphology is observed. This result solves a con¬ 
troversy concerning a natural initial bias towards prolateness 
shapes claimed by Heavens and Peacock (1985) and Dubinski 
(1992), or towards oblateness infered from the Bardeen et al. 
(1986) results. Moreover, the recent determinations on the in¬ 
trinsic shape of elliptical galaxies by Fasano & Vio (1991) and 
Lambas et al. (1992), which resulted consistent with triaxial 
objects, suggest that the distribution of shapes survives the 
non-linear evolution. 

We shall study the aspherical collapse of perturbations in 
the non-linear density field, investigating the relationship be¬ 
tween the initial asphericity and the stability of properties of 
the distribution of collapse structures in the presence of a tidal 
field and in isolation 
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